ISOMONODROMIC DEFORMATION THEORY AND THE 
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Abstract. In 1980 Jimbo and Miwa evaluated the diagonal two-point corre- 
lation function of the square lattice Ising model as a T-function of the sixth 
Painleve system by constructing an associated isomonodromic system within 
their theory of holonomic quantum fields. More recently an alternative isomon- 
odromy theory was constructed based on bi-orthogonal polynomials on the 
unit circle with regular semi-classical weights, for which the diagonal Ising 
correlations arise as the leading coefficient of the polynomials specialised ap- 
propriately. Here we demonstrate that the next-to-diagonal correlations of the 
anisotropic Ising model are evaluated as one of the elements of this isomon- 
odromic system or essentially as the Cauchy-Hilbert transform of one of the 
bi-orthogonal polynomials. 



For the square lattice Ising model on the infinite lattice an unpublished result 
of Onsager (see [13]) gives that the diagonal spin-spin correlation {aofiaN,N) has 
the Toeplitz determinant form 



(1) {(ToflCTN.N) = det(a^_j(A:))i<i.j<Ar, 

where the elements are given by 

d6 kcosnO — cos{n—l)6 



(2) 



— IT 



27r y/k^ + 1- 2fccos( 



A significant development occurred when Jimbo and Miwa (TUl [H] identified H]) as 
the T-function of a Pyi system. This identification has the consequence of allowing 
([T]) to be characterised in terms of a solution of the cr-form of the Painleve VI equa- 
tion, a second order second degree ordinary differential equation with respect to 
t := k^^ with parameter TV, or as the solution of coupled recurrence relations in TV 
with parameter t, which were subsequently shown to be equivalent to the discrete 
Painleve V equation. This was derived from the monodromy preserving deforma- 
tion of a certain linear system as a particular example of their general theory of 
holonomic quantum fields [14[ I17[ 112] , however the theoretical machinery employed 
there was never put to use on related problems arising from the Ising model. See 
the forthcoming monograph |16| on recent progress utilising this viewpoint. In a 
recent work [5] Forrester and the present author identified H]) as a r-function in 
the Okamoto theory of Pyi |15| and subsequently developed an alternative isomon- 
odromic theory jSj founded on bi-orthogonal systems on the unit circle with regular 
semi-classical weights. Wc remark that a result of Borodin [3] can also be used for 
the same purpose. 
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In a further development Au-Yang and Perk [2], [I] discovered that the next-to- 
diagonal spin-spin correlations have the bordered Toeplitz determinant form 



(3) 



(ccofw.Ar-i) = det 



ai 
\aN-i 



a-N+2 &Ar-l\ 



, N >1 



ai 



bo J 



where the elements a„ are the same as those above and the 6„ are given by 
(4) 

d0 C {kS - S)cosne + kScos{n-l)e - Scos{n+1)0 



27r Vfc2 + 1 - 2fccos6» 



S'2 + S'2 + 2fccos( 



(the definitions of the model parameters fc, S, S are given in the following para- 
graph). The task of the present study is to answer the following questions - can 
this correlation be evaluated in terms of a Painlevc-typc function and if so which 
one? The answer is in the affirmative and we identify the function in Proposition [3l 
In order to understand the result for the ncxt-to-diagonal correlations in its proper 
context we will need to revise some relevant known results for the diagonal correla- 
tions. In fact even in an algorithmic sense in order to compute the ncxt-to-diagonal 
correlations one has to first compute the diagonal ones. 

Consider the Ising model with spins (7r G {^1, 1} located at site r ~ (hj) on a 
square lattice of dimension (2L + 1) x (2L + 1), centred about the origin. The first 
co-ordinate of a site refers to the horizontal or x-direction and the second to the 
vertical or y-dircction, which is the convention opposite to that of McCoy and Wu 
|13j and early studies where the first co-ordinate labeled the rows in ascending order 
and the second the columns from left to right. We will focus on the homogeneous 
but anisotropic Ising model where the dimcnsionlcss nearest neighbour couplings 
are equal to K and K in the x and y directions respectively (see e.g. [3]). 



(*,.7 + l), 


1 1 


+ .7 + 1) 


K 


K 

1 ( 











Figure 1 . Co-ordinate system and couplings for the homogeneous 
anisotropic square lattice Ising model 
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The probability density function for configuration {o'ij}fj=^L given by 

^ L L-l L L-l 

(5) Pr({(Ty}^j=_i) = exp [i? a,ja,+ij+K ^ ^ (T^j(T,J+l , 

^^^^ j — — L i— — L i— — L j — — L 

and averages are defined by 

(6) (•)=E-P^({'^»^}^Wl)- 

The normahsation Z2L+1 is the partition function and conventionally periodic 
boundary conditions, cri,L+i = Ci^-L, o'l+ij = for all are taken for 

convenience. In all such averages the thermodynamic limit is taken limi_»oo(-) 
keeping K, K fixed. The relevant variables in our study arc the following variables 
k,S,S, (7,(7 defined by 

(7) 5* :=sinh2K, 5:=sinh2/?, C:=cosh2A', (5:= cosh 2^:, k := SS 

We will only treat the system in the ferromagnetic regime K,K > and k € 
(0,00), which exhibits a phase transition at the critical value k ~ 1. We will find 
subsequently that, from the point of view of the theory of isomonodromic systems, 
that the next-to-diagonal correlations are functions of the two complex variables, k 
and one of S, S, with k playing the role of the deformation variable and —S/S the 
spectral variable. While all of the results can be continued into the complex plane 
k, S € C suitably restricted we may often only state them for the physical regime 
k,S,S £ (0, 00). Corresponding to the Ising model is a dual partner Ising model, 
which is related to the original by the duality transformation or involution 

(8) CTr l-^ /^r, ("'ri • ■ • Cr„) 1-^ (/in • ■ • Mr„), 

(9) k^^, S^l, S^^. 

The dynamic variables fir are known as the disorder variables and can be given an 
interpretation in terms of the spins variables ar [12] . 

The appearance of Toeplitz determinants such as those of ([T]) is indicative of 
several structures and the most general of these is averages over the unitary group. 
Let U S U{N) have eigenvalues zi = e'^^ , . . . , zjy = e*^" . The unitary group U{N) 
with Haar (uniform) measure has eigenvalue probability density function 

^^^^ (2 )^Nl n \^k-zj\', zr.^e^'^eT, 9ie{-7r,7Ti 

^ ■ l<j<k<N 

with respect to Lebesgue measure dOi ■ ■ -(IOn (see e.g. [3 Chapter 2]) where T = 
{z € C : |z| = 1}. A well known identity [TH] relates averages of class functions, 
in particular products of a function w{z) over the eigenvalues, to the Toeplitz 
determinant 

N 

(11) I%[w] := (^Y[w{zi)z^'^ = det[w-,+j-kh,k=i,...,N, e € Z, TV > 1. 

1=1 

By convention we set /q = 1 and use the short-hand notation /jv ■= We identify 
w{z) as a weight function with the Fourier decomposition 

(12) w(z)=^u;iz'. 
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(14) a(C;fc) 



The specific Fourier coefficients appearing in the diagonal Ising correlations ([T]) 

are 

(^■^\ , f dC l l~k-^C-\ _ r d0 k cos n0^cos{n-l)e 

The implied weight is 

„^.-l/2^1/2(fc-l„^)-l/2(^„fc)l/2^ 0<fc<l ■ 

The analytic structure is different depending on fc > 1 (low temperature phase) or 
fc < 1 (high temperature phase). The reason for the phase change of e~'^^ in the 
weight is because of the argument changes 

(15) ^_fc = e-"(fc-C), fc-^-C = e"(C-fc-^), 

as fc goes from the fc > 1 to the fc < 1 regime. The correlation function for the 
disorder variables is 

(16) {^J■o,o^J'N,N) = det(aj_j(fc))i<ij<Ar, 
where 

^^7^ ~ - f ^" Z ^-^^'' _ T d0 cosn0 - kcosjn - 1)9 
^7) J^2m(:'' V 1-fc-C ^J-^2n Vk^ + l-2kcos0 ' 

The weight is 

^-fc■l/v/2(fc-c)-^/^(c-fc-l)^/^ i<fc<oo 

^ ' ^ " \fcl/2^1/2(C-fc■)"l/2(fc-l_C)l/^ 0<fc<l 

Although we use the same notation for the Toeplitz elements as Au-Yang and 
Perk [Tj the relationship between our elements and theirs is a„ — a^~^^ and 

From the viewpoint of the work [8] the weights (fT4|l . (fT8|) are particular exam- 
ples of the regular semi- classical class, characterised by a special structure of their 
logarithmic derivatives 

(19) 4T:f-W-wT = E^' p^^^- 

w{z) dz W{z) p^^z- Zj ■' 

Here V{z), W{z) are polynomials with degl^(z) < 3,degW^(z) = 3. The data for 
the weight (|14p is then 

(20) {z,}%, = {0, fc-\ fc}, {p,})^, = {V2, - V2, V2}. 
The data for the other weight is (flS)) 

(21) {z,}U = k-')^ {P^)U = {V2, -V2, V2}, 

which is the same as the previous case except for the exchange in the position 
of variable singularities. A particular observation in the Ising model case is that 
the Toeplitz matrix is not hermitian and the weight 'w{z) is complex for real and 
physical fc € (0, 00), z G T. The duality transformation is simply a transposition of 
the singular points Z2 ^ Z3 and at the critical temperature these two singularities 
coalesce. 
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An important identity relating the dual Toeplitz elements to the direct ones is 
the following well known duality relation [l] . 

Proposition 1. For all k and n we have 

(22) a„(fc) a„(fc~^) = -a_„+i(fc). 
The two weights are related by the duality transformation 

(23) a(C;fc)=a(C;fc-i). 

By regarding the Fourier integral in (jl3p as a contour integral, and changing the 
contour of integration one obtains the well known fact that the Toeplitz elements 
in the low temperature regime are given by 

(24) a„--^^^^i-/^J^^fc-"2i^i(V2,n-V2;n+l;fc-2), n > 0, 

(25) a_„ = ^^""^/^^l^/y^^ fc-"2j^i(-V2, n + V^; n + 1; k''), n > 0, 

7ri [n + 1) 

whilst those in the high temperature regime are 

(26) an ^ ' k^'-\F,{~%n~ y,;n;k% n > 1. 

TTL (n) 

(27) a^„ = '^1?^'^ fc"+^2f 1 ( V2, n + 1/2; n + 2; fc^) , n>-l, 

These elements are expressible as linear combinations of the complete first and 
second elliptic integrals K, E with arguments k~^ and k respectively [9] and with 
coefficients polynomial in these arguments. In the ensuing discussion we adopt the 
following shorthand notation for the complete elliptic integrals of the first kind 

(28) K< K(fc), K> K(fc-i), Ko := K(fco), 

with analogous notation for the second kind and where k<) = 2^/k/{k + 1) is the 
inverse Landen transformation. The complementary modulus is defined k'^ := 

The complex weight w(z) with support contained in T implicitly defines a system 
of bi-orthogonal polynomials {0„(z), 0„(z)}5^q on the unit circle by the orthogo- 
nality relation 

(29) I 4^MC,)<l^m{C,)4>n{C) = 5m,n, 



I J 2'iTiC^ 

whose existence is ensured if and only if /„ ^ for all n £ N. Notwithstanding the 
notation, (/>„ is not in general equal to the complex conjugate of 0„ and independent 
of it. The leading and trailing coefficients of these polynomials 

(f)n{z) = K„z" + . . . + (/>„(0), 
^„(z) = K„z" + . . . + (t>n{0), 

occupy an important role in the theory where again 0,i(O) are not in general equal to 
the corresponding complex conjugate. With the so-called reflection or Vcrblunsky 
coefficients specified by 

(30) r..:=m f„:=*ffi. 
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it is a well known result in the theory of Toeplitz determinants that 

(31) In+l[w]In-l[w] _ n^-Jl^ n>l 

Rather than dealing with it is advantageous to define the reciprocal polynomial 

^ni^) by 

(32) C(z):=z"0„(l/z). 

In addition to the polynomial pair </>„, </>* we require two non-polynomial solutions 
of the fundamental recurrence relations appearing in the theory [5], 

(33) e„(z) := f f.^^/ " wiQMO, n>l, 

(34) e:(z) 4^}-:^w{C)mi n>l. 



Kn Jt 27riC C - ^ 
These form a matrix system 



(35) y„(z;t) 



(f)niz) en{z)/w{z) 

Kiz) -e*(z)/u;(z) 



which, for regular semi-classical weights, has the property [5] that their monodromy 
data in the complex spectral z-plane is preserved under arbitrary deformations of 
the singularities Zj. 

From the Toeplitz determinant formula ([1]) we observe that 

(36) {cro.o<^N,N) = dct[aj_k]j.k=o....,N-i = ^w[a(C; k)] := Iwik), 

and apply the known results of Subsection 3.1 in [7] which provides the following 
recurrence scheme for the diagonal correlations. 

Corollary 1 ([7]). The diagonal correlation function for the Ising model valid in 
both the low and high temperature phases for N > 1 is determined by 

/„„s (o'0.0O'Ar+l.Af+l)(o-0,0O'7V-l,Af-l) , 

(37) = 1 - rNrN, 

along with the quasi-linear 2/1 

(38) (2iV -I- 3)(1 - rNfN)rN+i - 2N [k + k'^ + {2N - l)rjvfjv-i] rjv 

+ (2^ - 3) [(27V - l)rjvfjv + 1] rN-i = 0, 

and 1/2 recurrence relation 



(39) {2N + 1)(1 - rNfN)fN+i -2N[k + k'^ - {2N - 3)fArrjv_i] 

+ {2N - 1) [-(27V + l)rjvfjv + 1] f^-i = 0, 
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subject to initial conditions ro ^ tq ~ I and 



(40) 
(41) 
(42) 
(43) 




1 - fcp Eq 



T/ie initial values of the correlations are 



(44) 
(45) 



ao 



7rE>, 



1 < fc < oo 



^[(fc2-l)K<+E<], 0<fc<l 
! 1 



TT 1 



A consequence of the duality relation ([T]) are the following obvious relations 
amongst the coefficients of the bi-orthogonal polynomial system. 

Proposition 2. For all n and k we have 

(46) m = (-l)"/-^-na] 

(47) f„[, 



1 



rn a 



Now we turn our attention to the object of the present study - the evaluation of 
the next-to-diagonal correlations. Let us recall that the elements bn of the bordered 
Toeplitz determinant ([4]) can be written as 



(48) 



C 



lk/c-1 



lT27TiS + SC\ fcC-1 

These elements will also have complete elliptic function representations however for 
the anisotropic model we require the complete third elliptic integral defined by 



(49) 



n(?i, fc) 



1 



Vl-fc2 sin^^l-J^sin^ 



We also adopt a notational shorthand for these, analogous to that for the first and 
second integrals 

(50) n< := n(-s'2, fc"), n> := ^(-l/5^ fc^^), no := n(-4fc(5 - sy^, k(y). 

We note that IIo is not analytic at 5^ = 5* and in fact has a discontinuity there of 
the following form 



(51) 



-sgn[S - S) + 0[S - S), asS^S. 
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The first correlation in this sequence {N ~ 1) has the elUptic function evaluation 

^ [C^Uy - K>] , 1< A: < oo 
K<] , 0<A:<1 



(52) 



(CTonCTl.o) = ^0 = 




(53) 



2t:S 



{S^ - 1)K 



+ §6(5-5'), < fc < cx), 



5-5 

where Q{x) is the Heaviside step function. The term with the step function in ((53|) 
is necessary to compensate for the discontinuity in n<j as given in (|5ip in order 
that the correlation function remain continuous at 5 = 5. The second correlation 
function (N = 2) has the evaluation 



(54) 



AC 
PF5 



[fc2(l _ ^2)E^ ^ (^.2 _ l)52j^^] n> 



4T72 , n ;.2^e2T^2 , ,.2^c2 A;2)E>K>}, 



+ /fc4E2^ + (l-fc2)52K^+A:^(5' 
[(fc2 - 1)K< + (1 - 52)E<] n< 



l<k< 



El + (1 - fc2)K2< + (C252 - 2)E<K<}, 



C 1 + 



f [(1 - 52)Eo - fc^C^Ko] 



"■0 



5 + 5 
5-5' 



< fc < 1 

2tt 9(5-5) 



1 + fc;. 







CC 



(l + fc^)2 

The correlation functions for the disorder variables or dual correlations are given 



(55) 
by 

(56) (/io.oMAT.AT-i) = det 



—El + fc^(52 - S')Kl - (1 - 52)(1 - 52)EoKc 



where 
(57) 



CS 



\aN-i 

dC C 



d^N+2 ^N-lX 



iV > 1, 



n-1 



1 - fcc 



T27ri5 + 5CV l-fc/C 



The correlations in this sequence also have elliptic function evaluations analogous 
to (j52ll55p but we refrain from writing these down as they can be obtained from the 
direct correlations using the duality transformation 

(58) {l^O,Ol^N,N-l) — (o'0,0O'jV,JV-l) I k^l/k ■ 

S^l/S 
S^l/S 

In addition the (fTo^ofAf-i.Af) correlations can be obtained from ((To,of w.^v-i) under 
the exchange 5 5. 

These correlation functions are in fact characterised as a solution to an isomon- 
odromic deformation problem associated with the particular sixth Painleve system, 
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which itself characterises the diagonal correlation functions. This observation is the 
key result of the present study. 

Proposition 3. The next-to- diagonal correlation functions are given by the sec- 
ond type of associated functions appropriate to the weight ([I^ evaluated at a 
specific value of the spectral variable 

C In-1 



(59) 



-liz = -S/S), 



2S kn-1 ^ 

and valid for N > I. Here In and kjv are defined respectively by ill]) and i31]) 
appropriate to the weight |_?^[ ). 

Proof. A result in the general theory of bi-orthogonal polynomials is the detcrmi- 
nantal representation with a Toeplitz structure for the reciprocal polynomial [8] 

/ Wo ... w_n+i z"\ 



(60) 



Kiz) = -^det 



\ Wn ... Wi 1 J 

Using this and the definition of second associated function ([M)) one obtains an 
analogous bordered Toeplitz determinant [12] 

/ Wo 



(61) 

where 
(62) 



■det 



Wn 



-2z 



\ Wn 

dc 



W-n+1 gn\ 

wi go J 
■w{C), z^T. 



T 27r< ( - z 

The evaluation ([59]) then follows by comparison of these last two formulae with ([3|) 
and gHl). □ 

Many consequences flow from this identification - all of the general properties of 
the associated functions [S] can be applied. One particular useful characterisation of 
the next-to-diagonal correlations is that they satisfy a linear three-term recurrence 
relation. 

Corollary 2. The associated function I159\) satisfies the generic linear recurrence 
relation 



(63) 



Ul(^) 



r„+ize;_i(z) = [r„ -f- r„+iz]e; (z), 



subject to the two initial values for eQ,e| implied by 159\) and 15S153\54\55\) . The 
auxiliary quantities appearing in 163\) and I159\} satisfy the generic recurrences 

In 

9 ' 



(64) 
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subject to their initial values 

(65) lo = 1, 4^ —, 

ao 

utilising |^5[ j. 

Wc remark that this associated function also satisfies a hnear second order 
difi'erential equation in the spectral variable z whose coefficients are determined 
by the auxiliary quantities discussed above. However we refrain from writing this 
down as it doesn't appear to have as much practical ultility as the recurrences in 
the above Corollary. 

To close our study we examine a number of limiting cases, namely the zero 
temperature, the critical temperature and high temperature limits. At zero tem- 
perature, k —f oo, the solutions have leading order terms (N > 1) 

(~72)iV,-Ar - (72)^ , -JV / \ 1 

At the critical point, k = 1, we have a complete solution for the bi-orthogonal 
system. The polynomial coefficients have the evaluations 

«^7\ .2 _ r(Ar + 3/,)r(iv + i/,) _ 1 - _ , 

^^^> T^nTY) ' '^^""(2iV+l)(27V-l)' 

which is consistent with the well known result |13| 

(68) ^"°'"""'"> = nr(, + v,)r(,-V2)- 

The isomonodromic system is 

(69) 0.(.) = - (,^^;);^_,) • 2^i(%, -N; -N + 3/, 

(70) 0^(z) - ■ 2^1(1/2, -N; -N + 1/2; z), 

(71) -KNENiz) = — • 2^^1(3/2, N+1;N + 5/2; 1/^), 

\ I ^ N N\ ) (27V + 3)(2iV+ l)z ^ ' ^ /2, / 

(72) i^ive^(^) = 2^1(1/2, N + 1;N + 1/z). 

This last result ((72|) is consistent with the critical next-to-diagonal correlation given 

inia 

(73) {cTo,oaN.N-i) = (ao,oaAr,Ar)C • 2^1(1/2, iV;iV + 1/2;-^')- 
At infinite temperature, k 0, the leading order terms arc (A^ > 1) 

and the series expansion of these about A: = in terms of the generalised hypcrge- 
omctric function is given in [7]. 
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